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CONVERSION FACTORS, U.S. CUSTOMARY TO METRIC (ST) 
UNITS OF MEASUREMENT 


U.S. customary units of measurement used in this report can be converted 
to metric (SI) units as follows: 


Multiply by To obtain 
inches 25.4 millimeters 
2.94 centimeters 
square inches 6.452 square centimeters 
cubic inches 16. 39 cubic centimeters 
feet 30. 39 centimeters 
0.3048 meters 
square feet 0.0929 square meters 
cubic feet 0.0283 cubic meters 
yards 0.9144 meters 
Square yards 0.836 square meters 
cubic yards 0.7646 cubic meters 
miles 1.6093 kilometers 
square miles 259.0 hectares 
knots 1.8532 kilometers per hour 
acres 0.4047 hectares 
foot-pounds 1355 5,5 newton meters 
millibars LOO? 3 10-2 kilograms per square centimeter 
ounces O35 OS grams 
pounds 453.6 grams 
0.4536 kilograms 
ton, long 1.0160 metric tons 
ton, short 0.9072 metric tons 
degrees (angle) 0.1745 radians 
Fahrenheit degrees 5/9 Celsius degrees or Kelvins! 


OO _———————X—X—X—X—X—X—X—X—X—X—XX 
1To obtain Celsius (C) temperature readings from Fahrenheit (F) readings, 
USS Etoponmilee (Cs (5//9)) Gs S52). 


To obtain Kelvin (K) readings, use formula: K = (5/9) (F -32) + 273.15. 


EVALUATION OF THE COMPUTATION OF WAVE DIRECTION WITH THREE-GAGE ARRAYS 


by 
Dinorah C. Esteva 


I. INTRODUCTION 


Wave direction is an important parameter in the solution of many 
coastal engineering problems. A knowledge of wave direction is essen- 
tial for (a) estimating the direction and magnitude of sediment trans- 
port by waves, (b) using refraction calculations to infer wave conditions 
at one site from measurements made elsewhere, and (c) verifying theories 
of wave generation. 


Visual observations of wave directions have been collected by ship- 
board observers for over a century. About 20 years ago the Beach Ero- 
sion Board (BEB), predecessor to the Coastal Engineering Research Center 
(CERC), engaged the assistance of U.S. Coast Guard installations in the 
coliection of visual observations of breaker direction from shore. How- 
ever, objective determinations of wave direction are desirable without 
being restricted to location, time of day, or visibility condition. The 
capability to do so involves the use of wave measuring instruments. 
Panicker (1971, 1974) presents extensive reviews of reports dealing with 
the determination of wave direction from instrument records with partic- 
ular euphasis on those involving sea-surface elevation or pressure records. 


In March 1970, CERC installed an array of wave gages at Pt. Mugu, 
California. Records from the array were to be used to compare redundant 
values of wave direction and to estimate the level of accuracy of the 
computations. The available procedures for determining wave direction 
from an array involved assumptions that had not been thoroughly estab- 
lished. Thus, the records from the array would also be used in a syste- 
matic examination of these assumptions, and of the reliability of wave 
gages. 


This study discusses the array performance and the information gained 
about wave direction. Redundant values of directions were obtained from 
the 10 three-gage arrays possible with five gages. The mathematical model 
used assumes that the sea surface is the result of the superposition of 
a small number of narrow-banded wave trains consisting of long-crested 
waves traveling in well-defined directions. It was also assumed that 
only one wave train is present with a particular period. The first 
assumption is supported by the energy spectra computed at CERC (Thompson, 
1974), by aerial photos of the sea surface (Fig. 1), and by radar images 
of the wave field (Fig. 2). Many published reports include photos similar 
to that in Figure 1; e.g., McClenan and Harris (1975). Fujinawa (1974, 
1975) conjectured that narrow directional spread might be responsible 
for the incomplete recovery of the true directional spectrum from field 
records in his computations using high directional resolution. 


Average values of wave direction for bands of constant frequency 
width were computed from cross-spectra between gage pairs. Direction 
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Figure 1. Aerial photo of wave field at Pt. Saint George, California. 
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Figure 2. Radar scan of wave field at Nauset, Massachusetts. 


estimates for all bands 0.01 hertz wide between approximately 30 and 3 
seconds were obtained for the 10 arrays. The results displayed discrep- 
ancies of the order of 20° for those bands with central periods above 10 
seconds and of 180° for those with shorter central periods. It had been 
expected that the array would yield direction to better than 20° and for 
periods between 25 and 7 seconds. 


To isolate problems associated with the calculations, the propagation 
of narrow-banded wave trains across the array was simulated in a computer. 
The computational model was applied to the simulated observations using 
the maximum frequency resolution available from spectral computations 
based on 20-minute records. It was found that the directional results 
obtained with this model are highly dependent on the spectral width, 
both in frequency and direction, of the wave train involved and on the 
relationship between wavelength at the site and gage separations. The 
assigned directions were recovered within 1° for 16-second waves when the 
frequencies of the spectral components in the wave train differed by 0.003 
hertz or more, and the directions were spread within a 5° arc. This fre- 
quency separation results in a minimum difference in periods of 0.7 second 
for waves with periods near 16 seconds. 


Application of the same analyses to field wave pressure records with 
standard deviations above 0.61 meter (2 feet) resulted in an average dis- 
crepancy of 20° among computed directions for narrow-banded wave trains 
with periods longer than 10 seconds. Larger discrepancies resulted for 
shorter periods. Thus, accuracies no better than 20° can be expected for 
wave directions resulting from three-gage arrays. 


1. The System. 


A minimum of three gages is required for a unique determination of 
wave direction by most proposed models. Since these models make a few 
assumptions about the nature of ocean waves which have not been estab- 
lished, some redundance was thought to be necessary which would require 
a minimum of four gages. However, it was agreed that a five-gage array 
would increase the probability of redundance in the ocean environment. 
An array was designed at CERC by Leon E. Borgman, statistician-engineer, 
while on sabbatical leave from the University of California, Berkeley, 
when the experiment was being planned. He investigated the directional 
resolving power of several array geometries and concluded that the pat- 
tern shown in Figure 3 would be the most suitable for the conditions to 
be expected at Pt. Mugu (Borgman and Panicker, 1970). 


The array was installed off Pt. Mugu, approximately 80.47 kilometers 
(50 miles) northwest of Los Angeles (Fig. 4), in about 9.14 meters (30 
feet) of water, 0.76 meter (2.5 feet) from the bottom. The gages in the 
array are pressure transducers developed mostly at CERC (Williams, 1969). 
The heart of the system is a Fairchild pressure transducer which is potted 
inside a 2-inch Plexiglas tube (Fig. 5) (Peacock, 1974). A plastic tube 
filled with silicone oil transmits the pressure from the seawater to the 
pressure transducer. The silicone oil is separated from seawater by 
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rubber diaphragms. A Teflon scouring pad dipped in antifouling paint 
separates the rubber diaphragm from the end of the package which admits 
the seawater pressure. The instrument is mounted vertically on a tri- 
pod (Fig. 6). The signals from the pressure transducers are brought by 
cables to a recording and transmitting station onshore. 


Do Data Collection. 


The array went into operation on 27 March 1970. The water pressure 
at the five gages was registered continuously at a rate of four times a 
second on digital magnetic tape. Difficulties experienced during the 
first year with individual sensors were mostly of short duration and 
were presumed to be due to biological activity. However, major diffi- 
culties were experienced with the recording system, and on 16 March 1971 
the recorder was disconnected at Pt. Mugu and transferred to the CERC 
laboratory, then located in Washington, D.C. Recording from all five 
gages was reinitiated at CERC on 9 April 1971. Records were obtained 
continuously until 3 January 1972 when the recorder failed. During this 
period much of the data were useless because of an unacceptable level of 
high-frequency noise. The source of the noise was difficult to locate 
and was not eliminated from the signal until shortly before the failure 
of the recorder in January. Since 2 February 1972, records from three 
of the five gages at Pt. Mugu have been included in the time-shared 
recording of waves from east and gulf coast wave stations (Peacock, 
1974). In 1972, data were recorded for 20 minutes out of each hour; 
since February 1973, data have been recorded for 20 minutes out of each 
2-hour interval. 


II. FIELD DATA ANALYSIS 


The five gages in the array provided uninterrupted data for most of 
the first year of operation. Eight daily observations, each consisting 
of simultaneous 20-minute records from the five gages, were processed 
from these data. The observations had starting times within 1.5 hours 
of the weather synoptic times (0100, 0400, 0700, 1000, 1300, 1600, and 
1900 hours P.s.t.). The potential energy in the wave field is propor- 
tional to the variance of the time history of sea-surface elevation at 
a fixed location (Kinsman, 1965). For most conditions, the standard 
deviation of the surface displacement is one-fourth of the significant 
wave height. The standard deviation of the pressure at a fixed depth is 
roughly proportional to the wave height and may also serve as a measure 
of the wave height. 


The standard deviation of the recorded pressure was computed from the 
records of each of the five gages for eight records each day. The stand- 
ard deviation of the record from each gage was compared with the average 
of the five gages. If the standard deviation from any gage differed from 
the average by more than 20 percent, the record from that gage was deleted 
and the average recomputed from the remaining gages. A comparison of the 
individual standard deviations from the mean for that time period is shown 
in Table 1. 
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Table 1. Percent of observations where the largest 
departure of the standard deviations from 
the mean in the observations was as indi- 
cated (871 observations in 1970). 


Deviation from mean (pct) Percent of Observations 


41 
3 to 10 52 


11 to 20 


This comparison indicates the system operated consistently. The 
field wave records used for wave direction computation (discussed later) 
were chosen from the observations in Table 1, for which the standard 
deviations from all gages differed by 3 percent or less from their mean 


and for which the average significant wave height (uncompensated for 
attenuation with depth) was above 0.61 meter. 


Fourier analysis provides a reliable procedure for obtaining the 
periods of the most important waves. A Fourier analysis of sea-surface 
elevation (or pressure) with time results in the distribution of energy 
with frequency, usually referred to as the energy spectrum. The energy 
spectrum for the record from gage 5 in each of the eight daily observa- 
tions was computed using the Fast Fourier Transform algorithm developed 
by Cooley and Tukey (1967). The first 1,024 seconds of the 20-minute 


record was used in this computation. Gage 5 was chosen because of a 
good history of performance. 


Fast Fourier Transform computations yield the contribution to the 
variance at each of a set of frequencies which are harmonics of a funda- 
nental given by the inverse of the record duration, T. In this study, 
the frequencies of these harmonics are referred to as spectral frequenctes, 
and the corresponding periods, given by their inverse, as spectral pertods. 
“he energy content between 32 seconds and 3 seconds was used to normalize 
che spectrum. The lower limit on period was estimated from the thickness 
»f the water column above the pressure sensors. A summarized spectrum 
Fig. 7) was formed by combining the energy content in 11 adjacent spec- 
ral periods. The band width in the summarized spectrum was slightly 
arger than 10 * hertz (0.0107 hertz). The energy appearing at each 
pectral period in the pressure spectrum was compensated for attenuation 
ith depth by using the classical hydrodynamic pressure correction: 


_ cosh kh 
Bison) = cosh kAz ” (1) 


here k is the wave number, h the water depth, and Az 


the vertical 
istance of sensor from bottom. 


This resulted in a surface or compensated 
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Figure 7. Summarized pressure and surface spectra. 


spectrum. Surface root-mean-square (rms) values were obtained from the 
compensated energy at each spectral period. 


Since 1972, four simultaneous 20-minute records from each of the 
three gages included in the time-shared recording sequence are processed 
daily. More records are processed during special studies. The records 
processed from gages 1, 2, and 3 of the Pt. Mugu array start about 0100, 
0700, 1300, and 1900 hours (P.s.t.). The significant wave height, the 
distribution function, the first five moments of the distribution, and 
the pressure spectra are computed for these records in a study of wave 
record variability. These records are not analyzed for wave direction. 


1. Computation of Wave Direction. 


For a long-crested sinusoidal wave with frequency 6,7, propagating 
in direction a, (Fig. 8), the phase difference between locations 1 and 


2 with coordinates (x,,y,) and (xX5,y5) respectively, is given by: 
Ogg 2 AICS = &)) COS OF 2 OG = Yo) Sin 1S }) 5 (2) 

where k= 21/Le is the wave number associated with frequency 4,7, 

and L* is the wavelength. The subscript is used to indicate the 


possible presence of different wave trains with different frequencies 
and directions. 


The addition of the wave profile, n3, at a third noncolinear loca- 
tion allows solving for the sine and cosine of a. Thus, a unique solu- 
tion for the wave direction is obtained from the following equation when 
the signs of numerator and denominator are considered. (see App. A.): 


rien [(x, - 3) 95 - Gy - x) %,]/D a 
- [(v1 - Yo) 13 - (1 - 3) S12]/D] © 


where %,3 is the phase difference between the third and first locations, 
and D is a function of gage separation. 


Phase differences between locations for each different wave period 
are the only unknowns in the right-hand side of equation (3). Estimates 
of a representative phase difference between gage pairs for bands of 
constant frequency width are easily computed from cross-spectra of the 
wave (pressure) records. These spectra give, for each band, average 
values of the covariance of the wave records along two perpendicular 
directions. 


Substitution of these representative values of phase difference into 
equation (3) affords an expedient and economic means of obtaining esti- 
mates of a "representative" wave direction for each of the spectral 
bands, provided the results are of engineering use. The agreement among 
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Figure 8. Long-crested wave propagating in direction On. 


redundant computations of direction from the 10 three-gage arrays is an 
indication of the degree of confidence that can be placed on the result- 
ing directions. 


Directions were computed for a few observations (using the equation 
below) for the representative phase difference in each spectral band, 
where again, the sign of numerator and denominator must be considered. 


aK Quad, 
O12, = tan Conta (4) 


The subscript A in equation (4) indicates the value is for a given 
spectral band. The cospectrum, Co, and quadrature spectrum, Quad, 
are defined as: 


Con = F X1Xp COS $49 
Quad, = A X1X> SIN $)5 > (5) 


where x, and x, are the spectral amplitudes from wave records 1 and 
2, and R indicates summation over the adjacent spectral periods com- 


bined to make up a band. Jenkins and Watts (1968) showed that the defi- 
nitions in equation (5) are equivalent to the more standard definitions 
based on correlation functions. 


A computer output for runs using this approach is presented in Appen- 
dix B. Summaries of the results of the observations at 0700 hours, 25 
June, and 2010 hours, 28 June, are given in Table 2. The first column 
in the table gives the period at the center of the band; the second 
column gives the average percent energy in each band for the five gages. 
The last 10 columns give the computed "'representative'' direction of each 
band for the three-gage array. Directions are measured from the seaward 
normal with positive values counterclockwise and negative values clock- 
wise. The table shows that in these two observations disagreements in 
direction of the order of 20° are obtained for the longer period peaks 
and of 160° for the shorter period peaks. Examination of the computer 
output in Appendix B indicates the results are typical. 


The array had been expected to yield directions to better than XO 
for periods between 25 and 7 seconds. Understanding of the problems 
involved was sought by simulating the propagation through the array of 
narrow-banded wave trains traveling in a specified direction (discussed 
in next subsection). 


2. Simulated Data Analysis. 


In simulating the wave records for use, special consideration was 
given to wave period and to the difficulties arising in spectral analysis. 
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Table 2. Directional results with 0.01 hertz resolution for two observations. 


Period at |] Avg. pct 
center of 


band 


Three-gage arrays 


25 June 1970 at 0700 hr 


5299 44 | 156 50 42 46 145 | —131 | —147 |, —39 
6.40 32 | 157 | —44 | —44 | —44 147 135 | —153 30 
6.87 —23 | —45 | —33 | —43 | —39 | —149 | —136 | —166] —18 
7.42 56 | 150 35 50 42 150 24 46] —20 
8.06 33 | 162 3 | 155 | —49 8 5 —5 0) 
8.83 —12 —9 | -19 | -18 0 5 D; —1 0 
9.75 3 5 it tS) 10 12 16 18 
10.89 4 4 4 —2 16 —1 2} 8 ila 
12.34 it) 16 16 ih 27 14 14 20 Dil 
14.22 30 26 29 16 40 Sil 28 32 34 
16.79 25 22 Dall 10 33} 26 D5 Dail 28 
20.48 —1 0 0 —8 ol 26 25 Di, 28 
26.26 LUG |) Wad | Waly | dsl |} 120 178 144 109} 135 
28 June 1970 at 2010 hr 
5g) 7.44 169 52 | 136 | —38 |—31 | —28 |—148 | —136} —155)| —23 
6.40 10.85 —47 | —38 } —33 | —40 | —34 | —33 |—148 | —138 | —157} —25 
6.87 16.56 —40 | —34 | —29 | —36 | —30 | —27 |—-151] —22 46| —21 
7.A2 Oks) =A =39) | —=2.9) 5) 35130) |, — 27) | 152) = 20 48} —20 
8.06 19.61 =) |) 338) || =A) |) 315 |) 80) SA bya) 72 49| —21 
8.83 7.65 —37 | —29 | —24 | —32 | —27 | —20 —18| —15}| —14| —14 
O75 0.96 —16 | —15 | —12 | —16 | —16 —2 —3 —3 —4 0 
10.89 iLodL@) 0 0 0 1 —7 25 —3 8 al 24 
12.34 2.34 23 iL) 11 5) —4 26 14 14 14 14 
14.22 4.46 12 12 10 IZ 4 24 8 10 16 18 
16.79 1.61 19 14 1) 9 0 28 15 15 All Dall 
20.48 1.00 108 46 40 25 2 by 40 49 43 34 
26.26 0.70 —91 | —89 | —75 | —89 | —67 | —95 —92| -—66} —56| —72 


Assume a Sinusoid with frequency given as: 


2n(m + 6) 
TENNENT (6) 


where At is the interval of time between samples, |é| is less than 
or equal to 1/2, and NAt = 1, the record duration. 


Equation (6) provides for assigning frequencies which differ from the 
spectral frequencies. The contribution to the variance at the spectral 
frequencies of the sampled record is given by Sf as 


s2 = a2 + v2, (7 


where and by aces che} Lourterscoctincrents. 


aim 

Harris (1974) showed that for values of m near m (i.e., for spec- 
tral frequencies near the frequency of the sinusoid), and for m far 
removed from one and N/2, the approximations below are good estimates 
to the coefficients. 


: Asin 76 cos(> = 76) 


a > 
oe m(m - m+ 6) 
Pea 2 A sin 76 sin(o - 16) (8) 


m(m - m+ 6) 


Slow convergence of the energy toward the spectral period closest to 
the assigned period is clearly indicated. Thus, the energy is spread 
over adjacent spectral periods. This spreading, due to the finiteness 
of the record, is usually referred to as sptllover. 


The technique routinely used at CERC to decrease spillover is to 
apply the cosine bell data window as defined by: 


ie | 2mty \ 
WES 5 1 - cos T Wap BWI iso. oN 5 (9) 


where yz are the values in the original record. The Fourier coefficient 
for the resulting function yz, is given by: 


fe . A sin 16 cos(¢ - 16) 
ia S si p) 
2n(m - m+ 6)[(m- m+ 6)? - 1] 
AE A sin 76 sin(¢ - 76) 
by, = 


i 2n(m - m + 6)[(m - m+ 6)2 = 1] 


Thus, convergence is greatly increased and spillover is effectively 
reduced to three adjacent spectral periods. 
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The cosine bell data window was applied to the simulated records; 
therefore, it was sufficient to consider wave trains consisting of three 
sinusoids with nearby periods spread over at most six adjacent spectral 
periods. In general, the periods of waves in the ocean will differ from 
the spectral periods. Thus, the sinusoids in each simulated observation 
were specifically assigned periods differing slightly from the spectral 
periods by use of equation (6). Central periods of about 8 and 16 seconds 
were chosen to simulate 8- and 16-second swells. Swells with periods in 
this range are observed on the west coast. 


Each simulated wave record consisted of values of the wave profile 
at the five gage locations computed at 0.25-second intervals for 17.07 
minutes (1,024 seconds), to simulate the sampling rate and record dura- 
tion customarily used at CERC for field data. 


The three sinusoids were assigned specific directions and zero initial 
phase at the origin of coordinates, and were propagated across the array 
assuming a constant depth of 9.14 meters. Appendix C shows that the 
Finite Fourier Transform gives the correct phases for three sinusoids 
thus combined, provided the sinusoids are assigned the same direction, 
nearly the same amplitudes, and the frequency difference of each com- 
ponent to the nearest spectral frequency is the same. A frequency dif- 
ference of 0.134/1,024 = 0.000120 hertz was chosen. 


Characteristics of the first eight simulated observations are given 
in Table 3. Rough considerations of refraction using linear theory 
(McClenan, 1975) yield 22° from the normal to the coastline as the maxi- 
mum possible direction that waves with a 16-second period may have at 
the depth of the array. Directions within 21° N. and 21° S. from the 
normal were chosen for the first four simulated wave trains. Since waves 
with an 8-second period may approach the coastline at the array site from 
a much wider arc, directions up to 60° N. and 60° S. were used for the 
directions of the fifth to eighth simulated trains. 


The computed spectra for these eight simulated observations are shown 
in Appendix D (Figs. D-1 to D-8). In these figures the variance of the 
record, proportional to the energy, at each spectral period is plotted 
versus a linear frequency scale. No grouping of the variance at adjacent 
spectral periods has been made. These ungrouped spectra are referred to 
as high-resolution spectra. For spectra computed from 1,024-second records, 
this high resolution is approximately 0.001 hertz. 


The program used to compute these spectra is the same program used at 
CERC for the analysis of field data. Many spectra of field data computed 
with this program and summarized by grouping 11 adjacent spectral periods 
are given by Thompson (1974). 


The effect of spillover in the spectrum is shown in the figures of 
Appendix D. Each spectrum resulted from combining only three sinusoids; 
however, energy contributions appear at from five to nine adjacent spec- 
tral periods. 
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Table 3. Characteristics of simulated wave trains. 


Wave nd ee Bae 
| train (cm) 


11.43 All So 


! Approximate value. 
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The periods assigned to the sinusoids giving rise to simulated obser- 
vations 1, 2, 5, and 6, differed by exactly three spectral periods (see 
Table 3). The spectra for these four observations (Figs. D-1, D-2, D-5, 
and D-6) exhibit three maxima separated by two minima. The energy at 
these two minima may be interpreted as due to spillover since no sinu- 
soids were combined with the corresponding periods. 


The high-resolution spectra from field wave observations are dis- 
cussed later; however, the use of the minima in these spectra to esti- 
mate spillover and noise is discussed here. The average energy at the 
minima between 25 and 7 seconds in the high-resolution spectra of the 
field observations was used as a measure of spillover and noise. Only 
spectral periods displaying an energy content at least twice this "back- 
ground" energy were interpreted as possibly arising from physical wave 
components in the wave field. 


The average directions resulting from the 10 three-gage arrays are 
given in Table 4 (last colum). Only results of computations at the 
spectral periods closest to the assigned periods are shown. The table 
shows that for wave trains 1, 2, and 3 the computed directions for the 
10 arrays agree with the input directions to within 1°. 


The directional results for wave train 4 are correct only for spec- 
tral period 62. The main difference between this train and wave train 2 
is the narrower spectral width. Wave train 2 gave the correct directions 
for all spectral periods; wave train 4 did not. 


For simulated wave trains 5 to 8, the average directions seem meaning- 
less. To determine whether these poor results were due to a programing 
deficiency, another eight sets of simulated records were generated, inter- 
changing periods and directions. The computer output for the simulated 
observations is in Appendix D. This appendix and Table 5 should be 
referenced in the following discussion of additional simulated wave 
trains. 


The directional results for the sinusoids with periods clustered 
around 16 seconds were of the same quality regardless of the assigned 
direction. However, the directional results for the 8-second sinusoids 
indicate that the capability to sense the correct direction for these 
shorter waves depends on the orientation of the three-gage array rela- 
tive to the direction of propagation of the incoming wave. The result- 
ing directions, which differed by less than 87° from the assigned direc- 
tions, are given in Table 5. The top of each column in the table shows 
the shape and orientation of the array. These results are not surprising 
since the effective gage separation for the different gage pairs varies 
with orientation relative to direction of wave propagation. Table 5 
also shows that the more nearly equilateral arrays have wider direction 
discernability. The design considerations for the array indicated an 
effectiveness for wave periods between 25 and 7 seconds (Borgman and 
Panicker, 1970). 
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Table 4. Computational results at closest spectral frequencies for simulated wave trains. 


Wave Closest Spectral Amplitude Direction 
train spectral period 
frequency 
(1/1024 hz) (cm) (0) 
1 61 5.62 20S 
64 TD2 20S 
67 5.70 20S 
D 61 16.79 5.62 218 
64 16.00 Too) I5§ 
67 15.28 5.70 20 N 
3} 62 16.52 5) 5? 20S 
64 16.00 5.35 20S 
65 IS 37/5 1.19 20S 
4 62 16.52 baby 218 
64 16.00 5.35 258 
65 15.75 1.19 36N 
5 LAS 8.19 5.62 142N 
128 8.00 7.55 142 N 
131 7.82 5.70 142N 
6 D5 8.19 5.62 142N 
128 8.00 TESS 78N 
131 TES? 5.70 162 S 
U 126 8.13 5.52 142 N 
128 8.00 5.35 142 N 
129 7.94 1.19 142 N 
8 126 8.13 He) 142N 
128 8.00 7.14 76N 
129 7.94 5.83 I72.Se 
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The consistent and exact recovery of assigned directions achieved 
for the simulated 16-record wave trains is in part due to the use of 
the high computational resolution available in a large computer. Use 
of less exact data as available from recording instruments is expected 
to result in a less consistent and accurate recovery of the true direc- 
tion. Simulated observations 1 to 4 were rerun truncating the computed 
profile values to three digits as is commonly available from recording 
systems. The effect of this truncation is estimated to have introduced 
an error of the order of +0.127 centimeter (0.05 inch) in the instanta- 
neous values of the profiles. No appreciable differences in computed 
directions resulted by this truncation. 


3. Identification of Wave Trains from the High-Resolution Spectrum. 


A wave train in the real ocean is conceivably made up of several wave 
components with nearby periods propagating in approximately the same direc- 
tion. Simplistic idealizations of such wave trains are exemplified by 
Simulated observations 1, 3, 5, and 7. For wave trains 1 and 3, the 
long-crested wave model gave the correct direction of wave propagation 
for all 10 combinations (within 5°), not only at those spectral periods 
closest to the periods of the sinusoids combined, but at several adjacent 
ones on either side of these periods (see App. D, Figs. D-1 and D=-3). At 
some of these adjacent spectral periods, the contribution to the energy 
was several times the background level and nearly the same at the five 
gage locations. Thus, it can be assumed that a wave train in the ocean 
will give rise to a number of adjacent spectral periods in the high- 
resolution energy spectrum with energy content several times the back- 
ground level. This background level can be estimated by inspection of 
the minima in the energy spectrum, as discussed previously. A criterion 
for what energy level will be considered “high-energy content" can be 
set, and groups of adjacent spectral periods in the spectrum with high- 
energy content identified. These groups may each be assumed to arise 
from the presence of a wave train in the field with a mean wave period 
within the range of spectral periods in the group (a wave packet). The 
number of adjacent spectral periods in each group will be used as a 
measure of the width of the energy peak in the spectrum and indicates 
the spread in periods of the wave train. The spread in computed direc- 
tions at adjacent spectral periods in a group is an indication of the 
degree of directional organization in the wave train. Large spread in 
directional results may indicate the possibility that crossing wave 
trains with nearly the same period are present. As results for simulated 
wave train 4 indicate, the long-crested model based on the assumption of 
a single wave train at each frequency is not suitable for a determination 
of wave direction in such cases. Multiple wave trains at a single fre- 
quency may result from refraction around a shoal or island or from reflec- 
tion by a vertical wall. 


4. Spectra and Direction of Wave Propagation for Field Data. 


The energy and direction of wave propagation at each spectral period 
were computed for 44 field observations where the average uncompensated 
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significant wave height was over 2 feet and the discrepancy of individual 
standard deviations from their mean was 3 percent or less. Plots of the 
high-resolution spectra for these observations are in Appendix E. The 
vertical lines represent the energy contribution at each spectral period. 
The background level for each observation was estimated from the minima 
between 25 and 7 seconds in the spectra. Spectral periods in this range 
with energy content _above twice the estimated background energy were 
identified. Contributions to the energy satisfying this criterion at 
adjacent spectral periods were considered as arising from the same wave 
train. The number of adjacent periods in each train was used as a measure 
of the spectral width of the train. The energy had to be above the chosen 
level at all five gages for the spectral period to be included in the 
group. The spectral period among these showing maximum energy was taken 
as the "period" of the wave train. 


Directions were computed at all the spectral periods in each train 
for the 10 arrays. The total spread among these directions was found, and 
an average total spread was computed for the trains having the same spec- 
tral width. The same was done for the computed directional spread at the 
period of the train. 


Twenty-five percent of the identified wave trains had total direc- 
tional (computed) spreads above 100° and were not considered further. 
For 89 percent of the discarded trains, the period of the train was under 
9.4 seconds. Thus, all trains with periods under 9.4 seconds were dis- 
carded. 


Results for the different spectral widths for the trains retained 
(280) are shown in Table 6. The second column in the table gives the 
average total directional spread for the corresponding spectral width; 
the third column gives the average directional spread at the period of 
the wave train. The last column gives the number of wave trains having 
the spectral width in the first colum. 


These results indicate that the total directional spread increased 
with frequency width. Narrow peaks consisting of from one to three 
spectral periods are most frequent, and the spread in the direction at 
the period of the train remains relatively constant. Since the average 
spread for narrow-banded trains (width <0.003 hertz) is 21.8°, it is 
expected that three-gage arrays cannot yield directional results to any 
better accuracy. The mathematical exercise in Appendix C shows that 
array dimension is a limiting factor as to the shortest period for which 
some directional discrimination may be expected. An important factor in 
the validity of the directional result is the spectral structure of the 
wave train involved. Only in very special circumstances will the quanti- 
ties involved in equation C-10 (App. C) combine to give better results. 


There are various possible explanations for the large spreads observed 
in the directional results from field records. For the long-crested wave 
model to be strictly applicable, it is important that: 
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Table 6. Average spread in computed directions for 280 wave trains identified 
in the high-resolution spectra of 44 field wave observations. 


Avg. spread in Cases 
direction for 


Spectral width Avg. total spread 


in direction for 


period of train 


all periods 


(°) (No.) 
21.8 96 
22.9 58 
20.2 41 
20.4 22 
19.8 21 
19.8 16 
17.0 5 
17.6 8 
15.3 6 
ADA 3 
18.0 1 
15.0 1 
15.2 2 


(a) The phase differences be known accurately or that the 
probable error in their computed values be known. 


(b) The wave crests over the array site be long and 
straight; thus, the waves must not have undergone appreciable 
refraction. 


(c) The sea surface be stationary in time for the dura- 
tion of the record and in space over the span of the array. 


The mathematical exercise in Appendix D indicates that the analysis 
yields accurate phase differences only for strictly monochromatic condi- 
tions. When this is not the case, no accurate estimate of the error 
involved in the computation of direction can be given. This inability 
is inherent to the computational procedure and cannot be resolved. 


Waves with periods over 8 seconds have been and are undergoing 
refraction at the site of the array. Therefore, the wave crests are 
not exactly straight. For the longer waves, with wavelengths at the 
array site several times the gage separations, the curvature will not 
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introduce much error. This will not be the case for the shorter waves 
and orientation of the array becomes important. 


Because of refraction, the curvature of a wave train changes, perhaps 
only slightly, while propagating over the array. This change will intro- 
duce differences in the direction at each gage and differences in the 
direction sensed by different gage pairs, causing undetermined additional 
errors in the computation of direction. To determine the magnitude of 
these errors, two additional sets of simulated wave records were gene- 
rated. The periods of the sinusoids combined were those for simulated 
observation 3; the directions assigned were spread within a 10° arc for 
the first set and 20° for the second. The last two computer outputs in 
Appendix D show that spreads of the order of 16° and 32°, respectively, 
resulted in computed directions. 


A stationary condition in time is usually assumed when developing 
wave directional models. Indications are that this is not strictly 
applicable at all times. 


The three factors discussed above are sufficient to account for the 
inaccuracies encountered in the computations. 


5. Conclusions. 


The results of directional computations, for both simulated and field 
wave data records, indicate three-gage arrays have some capabilities to 
determine wave direction under certain conditions. These capabilities 
depend on: 


(a) The dimension of the array and the water depth at the 
Site which place a lower limit on the wave period for which 
possibly accurate directions may be computed. 


(b) The orientation of the array for the shorter periods. 


(c) The nature of the wave field; directional results for 
wave trains with a narrow frequency distribution or where the 
computed directions differ little at the adjacent spectral 
periods might be meaningful. 


For wave trains with narrow frequency and directional width, and 
period above 10 seconds, the three-gage arrays at Pt. Mugu yield direc- 
tions to an estimated accuracy of 20°. 


At the Pt. Mugu site, 16-second waves may approach the coastline at 


angles of 22° or less from the normal. The directional information pro- 
vided by the array adds little to this and seems hardly cost effective. 
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APPENDIX A 
DERIVATION OF THE EXPRESSION FOR WAVE DIRECTION 


Let the coordinates of nearby gage sites be Cepovgo 2S ip NE 
with N = number of gages. The water surface displacement at each site 
due to the passage of a sinusoidal wave of frequency o, and amplitude 
A, traveling in direction a (a7 in Fig. 8), is given by: 


pg SA cos {k [x, cos a + y, sin a] - 2not - $} , (A-1) 


where k = 27/L is the wave number, L the wavelength, and 4 the 
initial phase at the origin. 


The phase difference, ;;, between locations Iand sito the 
sinusoid considered is: 


ong S Iei(Osp x7) CS G > OR = yj) sina]. (A-2) 
Thus, for three noncolinear locations, 


ie S IOs > op) OS Ch 9 (67) 2 5%) Sab CH 


®13.= k[(x) = x3) cos o + (yy = ye) sin a], (A-3) 


suffice for a unique solution of the direction a. Eliminating first the 
sin a terms and then the cos a terms, to obtain: 


COs = 2s) Oia 2 es = 2a) Wale 


sin a 


k[ (x, ie X3) 4 = Vol) 2 (x; Xp) (YY a Y3)] 3 
Peston ch (Cra, Sia). Ong =, Oa 2 7g) One (A-4) 


kf Qa - x3) 71 - y2) - (1 - x2) (yi - y3)] | 


Since k is always positive, consideration need only be given to the 
other terms. Letting D stand for the quantity in square brackets, the 
direction). a) tor Dae Om ese canventsbiye 

on | LGn 2 26) 8na = Co = 2a) Ongl/o 


; A-5 
Gn = Ya) ©19 = On = sa) Gna] ee) 


a = tan 


A unique value for direction can be obtained by considering the signs of 
both numerator and denominator. The quantity D differs from zero for 
all nonlinear arrays as shown below. 


Let xX, = y; = 0, D will equal zero for yo/xX> = y3/x3; thus, 
Xo, Yo and x3, y,; will be on a straight line with slope given by 


this ratio. 
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APPENDIX B 


COMPUTER OUTPUT FOR CROSS-SPECTRA COMPUTATIONS 


Table. Guide to computer output. 
Explanation 


Title, plus date and time of the observation 
(day, month, year, hour, and minutes). 


Headings for columns. The numbers separated by 
dashes in the fourth line give the numbers of the 
gages in the array (see Fig. 3). 


Sequential number of bands 0.0107 hertz wide. 
Period at center of band (seconds). 

Percent of energy in each band for gages 1 to 5. 
Normalized to the energy content from approxi- 
mately 30 to 3 seconds. 

Resulting "representative" direction of wave 


propagation for the 10 arrays for the corre- 
sponding band (degrees). 
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APPENDIX C 


FOURIER COEFFICIENTS FOR A MIXTURE OF THREE SINUSOIDS 


The Fourier Transform of the function: 


f(t) = A cos(ot - 9) (C-1) 
with 
2n(m + § 1 
— : Is] sz > (cea) 


where 1 << m << N, computed from N values of f(t) evaluated at 
equal increments of t; At is given by the set of coefficients (Harris, 
1974): 


_ 2A sin 76 cos(mé6 - 6) 1 ie il 
a = N tan[n(m- m+ 6)/N] tan[n(m+ m+ 6)/N]]’ 
ne 2A sin 76 sin(md - 6) | 1 4 1 
is N tan[t(m - m+ 6)/N]  tan[t(m + m+ 6) /N] ; 


ma ty 2, 0 09 es (C-3) 


Harris shows that for values of m near m, and for m far removed 
from 1 and N/2, 


- A sin 76 cos(> - 76) 
a(m - m+ 6) 


a 


5 i 6 si - 16 
be 2 A sin 16 sin(¢ 16) aaa, a 


™(m - m+ 6) ? Cy 


Nz 


are good approximations to the coefficients in equation (C-3). Equation 
(C-4) shows that convergence is slow. 


For the special case 6 = 0, substitution into equations (C-3) or 
(C-4) gives: 
am = bm = 0, mf7m , 
and 


indeterminate for m=m. 


am = bm 
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Use of L' Hospital Rule in equation (C-3) when 6 +0, shows that: 
a, = A cos > , by, = A sin o 5 m=m. (C-5) 
Application of the cosine bell data window, 


£(nAt) == HE - cos a) f(nAt) ; Wie Ns 6 too Nog (C-6) 


is equivalent to replacing the original sinusoid f(nAt) with the sum 
of three sinusoids (Harris, 1974), where 


= [meet | 


m - 2n(m + 1+ 6 
- cog [G4 * Oe | o] - cos jpaaeee. ; en 


£(ndt) 


N| > 


For 6 = 0, and in view of equation (C-5), the Fourier Transform of this 
modified function will be given by: 


A [Nee 
any = ame = — 7 COs > > DA = Wang Tha oot 
and 
aj, = F cos Dba = 


with am = by, = 0 for all other values of m. Thus, energy appears at 
three adjacent m values. 


Harris (1974) shows that after application of the cosine bell data 
window, the approximate values of the coefficients are given by: 


A sin 16 cos(¢ -‘ 76) 


ia ane eR ERED TST RT PEGG ED 
2n(m - m+ 6)[(m- m+ 6)? - 1] 
é i i - 76 

op 2 aot Se ee ere Nes 
27(m - m + 6)[(m - m+ 6)* - 1] 2 


Thus , convergence increases rapidly, and values of the coefficients for 
(m - m) 2 3 may be disregarded. 


Equations (C-3), (C-4), and (C-8) imply that 


a) 


tan(@ - 16) = —”@ (C-9) 


m 


joy 
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Since 6 may be as large as 1/2, phase values computed from the 
Finite Fourier Transform may be in error by as much as 90°. 


For a simplistic simulation of a wave train, it is sufficient to com- 
bine three sinusoids with nearby periods propagating in the same direc- 
tion. Letting A; equal the amplitudes and k;, i- 1,2,3, the wave 
numbers, the Fourier Transform of the combination is given by: 


Az sin m6¢ cos(j - 157) 


3 
Doni. hp oick a EERE AGT eS LL oe 
Z, any o mo 8g) ly = me Sg)> oa 


i 
3 Af sin 167 sin(6@7 - 167 
beans L 4 (o7 “) : (C-10) 


a 


i=) 2m(mz - m+ 63) [ Gm; - m+ 63)? - 1] 
Nearby periods are attained by setting 
Nes (ie = mM) 6% % (C-11) 


Let the coordinates of three nearby locations be (xj,yj); j = 1,2,3. 
The only difference among the Fourier Transforms (eq. C-10) arising from 
wave records at each location is the values of the 6's. At each loca- 


av 
tion) ji, thes es) values are: 


on4 = ky (xj cos a + Yj sin a) - Oe i e=wl, 2.38 (C-12) 


or 
O77 = kiQz - oc , 
where 
Q3 = Xj cos a + yz sina , tubs Ap Sinn (C-13) 


Since the three sinusoids are assumed to have nearby periods, let 
eB iko 2 lke Sk 
Thus , 
One kQz - >4- 
Then, for the wave record at location j 
3 A; sin 76; cos(kQy - TO, — 7) 


amg = yp 9 
i=1 2n(d¢ + 62) [Ce + 62)? - 1) 


3 Ae sin 76” sin({kQ; - mo; - oO) 
eo weber Dees iE (C-14) 
fon Arle = OA) Wee Gao s iy 


90 


Let: 


au S Qin > Oe) Gs = Os)F ar ay) 


The expanded expression for ani after collecting terms in 
and sin KQ;, SPS 
: Ay Sin 16] cos(¢1 + 16j) 
ang = COs ko; i] 
Ag sin 152 coS(¢2 + 162) 
[2] 
A3 Sin 763 cos($3 + 163) 
i [3] | 
A; Sin 75, sin(o, + 176,) 
cp Sali kQ5 a 
Ay Sin 75 Sin(¢5 + 165) 
i [2] 
A3 sin 763 Ssin($3 + 763) ‘ 
: [3] 
and similarly for Ding 
A, sin 75, cos($, + 16)) 
Dig = sin kQ; | [a] 
Ay Sin 755 cos(o5 + 155) 
i [2] 
A3 Sin 753 cos($3 + 143) 
i [3] 
A; Sin 76, sin($, + 16)) 
- cos k{Qj | [1] 
Ap Sin 155 sin($5 + 155) 
: [2] 
A3 Sin 763 sin(¢3 + 153) 
qP [3] ° 


For every narrow-banded wave train, JA, | < 3 in equation (C-15), 
unpredictable terms are introduced in the expressions for the coefficients 


5 | 


which make the ratio a poor estimator of the value of the phase 


of the sinusoid. 


Din/ Ayn 


The numerators of the terms inside the braces in the above equations 
for ang and bing Boge BNE TOS Oi Omlor Mes a S hz 5ss 


Randomness in the values of 
tionships of the three sinusoids 
duce partial cancellations among 
the resulting error. 


é6;, i = 1,2,3 and in the phase rela- 
at the origin of coordinates might pro- 
the terms inside the braces to reduce 


Assume for example: $; = 0, 1 = 15255. The coefficients reduce to: 
aes A, sin 276 A> sin 276 Az sin 216 
ogg ~ a CE kQ; | zs Sy A a ae he ge mes 
[1] [2] [3] 
A, sin? 16 he Sim in2 
Bite ka; 1 75 4 2 Sin“ 15 i A3 sin* 163 
[1] [2] [3] 


Similarly, for Ding 


oi Gags A, sin 276, A» sin 2165 A3 Sin 2163 
Da ay SUR a re 
2 J [1] [2] [3] 
A, sin? 16 A, sin? 16 A, sin? 16 
- cos kQ; jp ee g S ae (C-16) 
[1] [2] [3] 
Letting: 
Ly = Ay sin 275 1 ’ Lo = Ay sin? 70 4 > 
M, = Ap sin 276, , M, = Ay sin* 6, , 
Ny = A3 sin 2163 , Ny = Az sin 13 ; (C-17) 
then: 
oO 2 1 A Ly My Ny 5 E Lo M N 
oe soz ae Ee “oT ml im om 
and 
6 ; L M N L 
ee ecarin el Ny OU. ese te (2), eae) 
2 ful wl fs] Lo) aS] 
Assume further: A, ~ Ay ~ A3 and | in, =ml o@ woe A= 2 anal equal 
for i= land 35. Sincem S| <5) the terms) [ll sande: |S) seanclor 


a2 


approximately the same magnitude but of opposite signs. 
involving the products [1] [2] and [2] [3] tend to 
6, = 65 = 63=6: 


Thus, terms 
cancel. Letting 


al sin kQ; sin 274 - 2 cos KQ; sin? 16 
eet = Wae IO. Gin Dae = > aim ikOo Gia2 om S 
ang cos kQ; sin 276 + 2 sin kQ; sin2 176 


Using the trigonometric identities for the double arc, 
reduces to: 


b . 
—WY = tan (kQ- - 18) 
ane J 


this expression 


(C-19) 


Phase differences between locations will be approximately correct. 


oe) 


APPENDIX D 


SPECTRA PLOTS AND COMPUTER OUTPUT FOR SIMULATED OBSERVATIONS 


Figures D-1 to D-8 show high-resolution spectra 
pressure gages 1 to 5 at Pt. Mugu, California. 
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Energy Density (cm?/ Hz) 
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Figure D-1. 
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Energy Density (em*/ Hz) 
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Figure D-2. 
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Energy Density (cm?/ Hz) 
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Figure D-3. 
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APPENDIX E 


HIGH-RESOLUTION SPECTRA FOR FIELD WAVE DATA 


Figures E-1 to E-44 show high-resolution spectra for 
pressure gages 1 to 5 at Pt. Mugu, California. Date 


and significant wave height are indicated for each 
figure. 
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Figure E-1. Significant wave height 93.6 centimeters (3.1 feet), 
10 April 1970, 2206 hours. 
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Figure E-2. Significant wave height 84.3 centimeters (2.8 feet), 
20 April 1970, 0928 hours. 
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Figure E-3. Significant wave height 81.0 centimeters (2.7 feet), 
20 April 1970, 1228 hours. 
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Figure E-4,. Significant wave height 105.3 centimeters (3.5 feet), 
21 April 1970, 0028 hours. 
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Figure E-5. Significant wave height 95.2 centimeters (3.1 feet), 
21 April 1970, 0229 hours. 
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Figure E-6. Significant wave height 100.7 centimeters (3.4 feet), 
21 April 1970, 0929 hours. 
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Figure E-7. Significant wave height 94.7 centimeters (3.1 feet), 
21 April 1970, 1229 hours. 
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Figure E-8. Significant wave height 113.8 centimeters (3.7 feet), 
21 April 1970, 1829 hours. 
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Significant wave height 85.1 centimeters (2.8 feet), 
21 April 1970, 2129 hours. 


88 


Period (s) 
25.0 15.0 10.0 6.0 


Energy Density (cm?/ Hz) 


0.06 0.10 0.14 0.18 
Frequency ( Hz) 


Figure E-10. Significant wave height 109.3 centimeters (3.6 feet), 
10 June 1970, 1421 hours. 
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Figure E-11. Significant wave height 117.9 centimeters (3.9 feet) , 
10 June 1970, 1731 hours. 
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Figure E-12. Significant wave height 113.1 centimeters (3.7 feet), 
10 June 1970, 2041 hours. 
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Falgune E13.) Slemita cant) wave hengh tus, jomcencimeterss (29 s-keet) i, 
10 June 1970, 2351 hours. 
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Figure E-14, Significant wave height 93.3 centimeters (3.1 feet), 
11 June 1970, 0301 hours. 
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Figure E-15. Significant wave height 81.0 centimeters (2.7 feet), 
11 June 1970, 0631 hours. 
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Figure E-16. Significant wave height 88.6 centimeters (2.9 feet), 
11 June 1970, 0942 hours. 
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Figure E-17. Significant wave height 71.0 centimeters (2.3 feet), 
24 June 1970, 0906 hours. 
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Figure E-18. Significant wave height 69.3 centimeters (2.3 feet), 
24 June 1970, 1206 hours. 
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Figure E-19. Significant wave height 79.5 centimeters (2.6 feet), 
24 June 1970, 1506 hours. 
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Significant wave height 72.5 centimeters (2.4 feet), 
24 June 1970, 1807 hours. 
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Figure E-21. Significant wave height 70.1 centimeters (2.3 feet), 
25 June 1970, 0007 hours. 
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Figure E-22. Significant wave height 76.9 centimeters (2.5 feet), 
25 June 1970, 0307 hours. 
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Figure E-23. Significant wave height 66.7 centimeters (2.2 feet), 
25 June 1970, 1808 hours. 
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Figure E-24. Significant wave height 70.3 centimeters (2.3 feet), 
25 June 1970, 2108 hours. 
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Figure E-25. Significant wave height 61.2 centimeters (2.0 feet), 
26 June 1970, 1758 hours. 


Period (s) 
25.0 15.0 10.0 6.0 


| 


Energy Density (cm?/ Hz ) 


| | | 


1 
26 0.10 0.14 0.18 
Frequency (Hz) 


Figure E-26. Significant wave height 79.7 centimeters (2.6 feet), 
26 June 1970, 2358 hours. 
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Significant wave height 83.2 centimeters (2.7 feet), 
28 June 1970, 2310 hours. 
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Figure E-29. Significant wave height 78.7 centimeters (2.6 feet), 
29 June 1970, 0510 hours. 
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Figure E-30. Significant wave height 64.0 centimeters (2.1 feet), 
16 November 1970, 1607 hours. 
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Figure E-31. Significant wave height 80.9 centimeters (2.7 feet), 
16 November 1970, 1907 hours. 
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Figure E-32. Significant wave height 63.0 centimeters (2.1 feet), 
17 November 1970, 0407 hours. 
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Figure E-33. Significant wave height 74.6 centimeters (2.4 feet), 
16 December 1970, 1700 hours. 
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Figure E-34. Significant wave height 95.3 centimeters (3.1 feet), 
16 December 1970, 2000 hours. 
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Figure E-35, Significant wave height 85.3 centimeters (2.8 feet), 
16 December 1970, 2300 hours. 
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Figure E-36. Significant wave height 90.9 centimeters (3.0 feet), 
17 December 1970, 0200 hours. 
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Figure E-37. Significant wave height 103.4 centimeters €3.4 feet), 
17 December 1970, 0500 hours. 
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Figure E-38. Significant wave height 120.7 centimeters (4.0 feet), 
17 December 1970, 0800 hours. 
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Significant wave height 97.4 centimeters (3.2 feet), 
17 December 1970, 1044 hours. 
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Figure E-41. Significant wave height 91.2 centimeters (3.0 feet), 
17 December 1970, 1644 hours. 
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Figure E-42. Significant wave height 108.1 centimeters (3.5 feet), 
17 December 1970, 1944 hours. 
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Significant wave height 71.5 centimeters (2.3 feet), 
18 December 1970, 0445 hours. 
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Figure E-44. Significant wave height 66.9 centimeters (2.2 feet), 
18 December 1970, 0745 hours. 
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